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ABSTRACT 

We examine the X-ray spectra of VII Zw 403, a nearby low-metallicity blue compact dwarf 
(BCD) galaxy. The galaxy has been observed to contain an X-ray source, likely a high mass 
X-ray binary (HMXB), with a luminosity of 1.3 — 23 x 10^® erg s“^ in the 0.3 — 8 keV energy 
range. A new Suzaku observation shows a transition to a luminosity of 1.7 x 10“^° erg s“^ 
[0.3 — 8 keV], higher by a factor of 7 — 130. The spectra from the high flux state are hard, 
best described by a disk plus Comptonization model, and exhibit curvature at energies above 
5 keV. This is consistent with many high-quality ultraluminous X-ray source spectra which 
have been interpreted as stellar mass black holes (StMBH) accreting at super-Eddington rates. 
However, this lies in contrast to another HMXB in a low-metallicity BCD, I Zw 18, that 
exhibits a soft spectrum at high flux, similar to Galactic black hole binaries and has been 
interpreted as a possible intermediate mass black hole. Determining the spectral properties 
of HMXBs in BCDs has important implications for models of the Epoch of Reionization. 
It is thought that the main component of X-ray heating in the early universe was dominated 
by HMXBs within the first galaxies. Early galaxies were small, metal-deficient, star forming 
galaxies with large H I mass fractions — properties shared by local BCDs we see today. 
Understanding the spectral evolution of HMXBs in early universe analogue galaxies, such as 
BCDs, is an important step in estimating their contribution to the heating of the intergalactic 
medium during the Epoch of Reionization. The strong contrast between the properties of the 
only two spectroscopically studied HMXBs within BCDs motivates further study on larger 
samples of HMXBs in low metallicity environments in order to properly estimate the X-ray 
heating in the early universe. 
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1 INTRODUCTION 

The first X-ray sources in the early universe were binary sys¬ 
tems containing compact objects formed in the collapse of mas¬ 
sive stars in low metallicity environments. During the Epoch 
of Reionization, X-rays from these high mass X-ray binaries 
(HMXBs) may have contributed to the heating of the intergalac¬ 
tic medium (IGM) (Shull & van Steenberg 1985; Haardt & Madau 
1996; Mirabel et al. 2011; McQuinn 2012; Mesinger et al. 2013). 

Recent models (e.g., Venkatesan et al. 2001; 
Ricotti & Ostriker 2004; Linden et al. 2010; Fragos et al. 2013) 
and observations (e.g., Mapelli et al. 2010; Kaaret et al. 2011; 
Prestwich et al. 2013; Basu-Zych et al. 2013; Brorby et al. 2014) 
suggest that X-ray emission and binary populations are enhanced 
(for a given star formation rate) in regions of low-metallicity. 
The theoretical studies find that: (1) lower metallicity pro¬ 
duces weaker stellar winds, which result in a greater number of 
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more massive black holes, and thus more massive black hole 
HMXBs (Fragos et al. 2013); and (2) that the HMXB popula¬ 
tions in low metallicity environments are dominated by Roche 
lobe overflow HMXBs which are very effective in transferring 
mass to the compact object in comparison to the wind-driven, 
supergiant HMXBs that dominate in high metallicity environ¬ 
ments (Linden et al. 2010). The observational studies find that the 
X-ray luminosity function describing the population of HMXBs 
remains relatively flat (dN/dL oc L~°‘, where a = 1.6), and 
therefore the total X-ray luminosity of the star-forming region 
is dominated by the most luminous sources. The bright end 
of the X-ray luminosity function consists of the aptly named 
ultraluminous X-ray sources (ULXs). 

Ultraluminous X-ray sources are bright (Lx > 10^® erg s-i) 
extragalactic sources that are not associated with the nucleus of 
their host galaxy. Studies have concluded that ULXs are likely bi¬ 
nary systems containing either a neutron star (Bachetti et al. 2014), 
an intermediate mass black hole (IMBH) (Colbert & Mushotzky 
1999; Kaaret et al. 2001; Dewangan et al. 2006), or a stellar 
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mass black hole (StMBH) with a luminosity that exceeds its 
Eddington limit (Poutanen et al. 2007; Gladstone et al. 2009; 
Kawashima et al. 2012; Sutton et al. 2013). The IMBH ULXs 
would have spectra similar to the StMBH binary systems within our 
Galaxy, but with higher luminosities (Colbert & Mushotzky 1999). 
However, if the compact object accretes at super-Eddington rates, 
then its spectrum will likely differ from the standard sub-Eddington 
black hole binary (BHB) states. Gladstone et al. (2009) provide an 
observationally defined state for a set of high quality ULX spectra 
called the ultraluminous state. This state is characterized by a soft 
excess modelled by a cool disc with a power law tail and a break in 
the spectrum above ~ 3 keV. The curvature at higher energies is not 
seen in the standard states of Galactic BHBs and Gladstone et al. 
(2009) identify this state with super-Eddington accretion flows onto 
StMBHs. 

The shape of the emission spectrum of the first X-ray sources 
plays a key role in understanding their contribution to heating of the 
IGM in the early universe and modifies constraints on HMXB pop¬ 
ulation in the first galaxies. The spectral shape affects how X-rays 
heat the IGM because the mean free path and the fraction of energy 
deposited as heat versus ionization varies depending on the X-ray’s 
energy (McQuinn 2012; Mesinger et al. 2013; Pacucci et al. 2014). 
The distance the X-rays travel and the energy they deposit in the 
IGM would affect the morphology of the heating and ionization 
during the Epoch of Reionization and these effects may be seen 
in future observations of the redshifted 21-cm hydrogen line (e.g., 
SKA^). 

Another effect of spectral shape is the modification of the K- 
correction for high redshift objects. If the most luminous X-ray 
sources in the first galaxies exhibit the same type of spectral curva¬ 
ture seen in ultraluminous X-ray sources (ULXs), as opposed to a 
straight power law, relatively fewer high energy photons would be 
redshifted into the soft X-ray band. This would relax the constraint 
imposed by the observed soft X-ray background allowing for a sig¬ 
nificant enhancement in the X-ray luminosity of early, star-forming 
galaxies (Kaaret 2014). 

In order to determine the spectral shape of ULXs in the early 
universe, we must study nearby analogues of the first galaxies. Star- 
burst galaxies, such as blue compact dwarf galaxies (BCDs), have 
X-ray fluxes that are dominated by the most luminous high-mass X- 
ray binaries within them (Thuan et al. 2004). BCDs are physically 
small galaxies with blue optical colors, low-metallicities, large H I 
mass fractions, and generally exhibit a recent burst of star forma¬ 
tion. These dwarf galaxies provide the best available local labo¬ 
ratories for young galaxies in the early universe (Kunth & Ostlin 
2000 ). 

Linden et al. (2010) showed that for Z/Z© < 0.2 the HMXB 
population of a star forming region should become dominated by 
Roche lobe overflow HMXBs. By extension, this would lead to a 
larger ULX population. Prestwich et al. (2013) used observations 
to show that the ULX population is enhanced at low Z/Z© and 
the transition takes place in the range 0.06 < Z/Z© < 0.13 
(7.7 < 12 -f log(0/H) < 8.0). The galaxies that were below this 
range (Z/Z© < 0.06) were deemed extremely metal poor galaxies 
(XMPGs). It is this population of XMPGs that we chose to focus 
on, as the earliest galaxies should be very metal poor. 

To date, a total of eight low-metallicity BCDs have been found 
that contain ULXs (Thuan et al. 2014) as defined by having a lumi¬ 
nosity > 10®® erg s“^ (Eddington luminosity of a 10 M© black 
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hole). Of these, six are XMPGs and only two of these BCDs, 
VII Zw 403 (Z/Z© = 0.019) and I Zw 18 (Z/Z© = 0.062), 
have had sufficient counts to allow for a spectral study. Other 
ULXs have been observed in sub-solar metallicity dwarf galax¬ 
ies (e.g., Holmberg II: Zezas et al. 1999; Kaaret et al. 2004; Holm- 
berg IX; La Parola et al. 2001; Wang 2002; Miller etal. 2004), 
and have metallicities (Z/Z© = 0.1 — 0.3 (Egorov et al. 2013) 
for Holmberg II and Z/Z© « 0.09 — 0.5 (Moustakas et al. 
2010; Makarova et al. 2002) for Holmberg IX) that may fall 
within the range for enhanced ULX production (Linden et al. 2010; 
Prestwich et al. 2013). However, for this study we strictly examine 
only the lowest metallicity galaxies, VII Zw 403 and I Zw 18, as 
these are the best analogues to the extremely metal poor galaxies in 
the early universe. 

In this paper we focus on VII Zw 403, a low-metallicity 
BCD that contains a luminous X-ray binary and is at a distance 
of 4.34 Mpc (Karachentsev et al. 2003). A recent Suzaku observa¬ 
tion of VII Zw 403 shows the X-ray source to be in a much higher 
flux state than previous observations and the spectrum is similar to 
many ULX spectra. In Section 2, we describe the procedures used 
to prepare and analyze our dataset. Results of timing and spectral 
analysis are presented as well. We explain the physical motivation 
for some of the more complex models and discuss them in the con¬ 
text of previous high-quality spectral studies. Within Section 3, we 
contrast the spectral evolution of VII Zw 403 with that of I Zw 18, 
another low-metallicity BCD containing an HMXB that has been 
observed to transition between low and high flux states. The ob¬ 
served spectral state transitions for I Zw 18 are more analogous 
to the Galactic BHBs than the ULX spectra, which is inconsistent 
with the observed ULX-like spectra of VII Zw 403. Larger spectral 
studies and more frequent observations of ULXs within XMPGs 
will allow for a clearer picture to emerge regarding the spectral 
states exhibited by HMXBs in the early universe. Throughout the 
paper all fluxes and luminosities are in the 0.3-8 keV energy range, 
unless otherwise specificied. 

2 OBSERVATIONS & ANALYSIS 

2.1 Timeline of VII Zw 403 X-ray Observations 

VII Zw 403 was observed with Chandra on 2000-01-07 for a 
total of 10.6 ks and about 300 counts were detected. The data 
were analyzed by Ott et al. (2005) who found a good fit with 
a hard power law model (T = 1.75) and an absorbed flux of 
1.5 X 10“®® erg cm“® s“® in the 0.3-8.0 keV band. This flux is 
lower by a factor of 40 than the high, hard state seen in the Suzaku 
observation discussed in this paper. 

There have been a total of 8 observations with Swift, two 
in 2011 and six in 2014. Data from 2011-01-27 and 2011-02-01 
show no significant source counts and are limited to fluxes be¬ 
low 4 X 10“®^ cgs. Table 1 and Eigure 1 show that the data 
taken from 2014-03-16 to 2014-04-06, roughly 4-5 months after 
the Suzaku observation, have flux levels varying around an average 
of 6 X 10“®® cgs which is 4 times higher than the low state (Chan¬ 
dra) and 10 times lower than the high state (Suzaku). There were 
not enough counts in the Swift observations to extract individual or 
summed spectra. 

VII Zw 403 was previously observed with XMM-Newton on 
2004-10-07 for 46.9 ks. We created new events files using the 
epproc and emproc commands with the most recent calibra¬ 
tion files. Events were selected from intervals with low background 
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Figure 1. Flux measurements of VII Zw 403 in units of erg cm~^ s“^. The 
triangles (green) are Swift data, star (black) for Chandra, circle (red) for 
XMM, and square (blue) for Suzaku. This plot is a visual representation of 
the data in Table 1. The inset shows the six Swift observations from 2014. 


Table 1. Swift observations of VII ZW 403. Count rates and absorbed fluxes 
are for the 0.3-8 keV band. Fluxes from Chandra, XMM, and Suzaku are 
shown for comparison. 


which (XISl) is a back-illuminated (Bl) CCD chip and the other 
two (X1S0,3) are front-illuminated (FI). The observation was made 
on 2013-12-01 with a total exposure time of 88.66 ks and resulted 
in 68409 net counts. 

The data were reprocessed using the standard HEASARC 
script aepipeline. Spectra for each of the three detectors 
were extracted from a circular region centered at RA(2000) = 
Ifh 2 gm Q^sg^ ^2000) = +78° 59' 37'.'1 with a radius of 
3' (90 per cent encircled energy fraction) and a background region 
from an annulus with inner and outer radii of 4' and 6', respectively. 
The spectral data were grouped such that there was a minimum of 
40 counts per energy bin. 

The spectra were fitted using multiple models in the 0.4 — 
10 keV energy range. For sufficient sensitivity, we fit the XlSl (Bl) 
detector data in the 0.4 — 8.0 keV energy range and the X1S0,3 data 
in the 0.6 — 10.0 keV range. 

Contamination has been building up on the spacecraft, leading 
to additional absorption, mainly by carbon. This effect can be seen 
in comparing identical spectra using two different sets of response 
files. Response files based on older calibration data (2006) would 
show more absorption in the models (wabs, tbabs, phabs, etc.) than 
response files generated with this contamination taken into con¬ 
sideration. For this reason, we generated our own response files 
for each spectra as opposed to using archived response files from 
the HEASARC archive. The response files were created by running 
xisrmfgen and xissimarfgen upon extraction of the spectra. 


Date 

Exposure 

(s) 

Net count rate 
(10-3s-i) 

Flux 

(10“^^ cgs) 

2011-01-27 

4207 

< 0.7 

< 0.4 

2011-02-01 

1664 

< 0.6 

< 0.3 

2014-03-16 

2851 

12.6±2.5 

6.3 

2014-03-22 

850 

20.0±6.1 

10.0 

2014-03-23 

753 

5.3±4.2 

2.7 

2014-03-28 

801 

12.5±5.3 

6.3 

2014-04-01 

1450 

9.7±3.3 

4.8 

2014-04-06 

1877 

12.3±3.1 

6.1 

2000-01-07 

Chandra 


1.5 

2004-10-07 

XMM 


0.5 

2013-12-01 

Suzaku 


60.0 


(±3cr of the mean quiescent level). From this, we found GTls of 
14.9 ks, 33.7 ks, and 32.3 ks for the PN, MOSl, and MOS2 de¬ 
tectors, respectively. Energy spectra were extracted from a circular 
region with a radius of 20 arcseconds and had a total of 630 net 
counts from the three detectors. The spectra were rebinned such 
that each energy bin had at least 15 net counts or was 1/4 of the local 
FWHM. The fitting was restricted to the 0.3-5.0 keV energy range 
due to background contamination above 5 keV. We fitted the spec¬ 
tra with a fixed Galactic absorption of Nh ~ 3.91 x 10^° cm~^ 
and a power law. The resulting fit led to a x^/D.o.F. = 35.5/33 
with a photon index ofF = 2.1 ±0.2 and a model flux of 
4.6 X lO”^'^ ergs cm“^ s“^ (0.3-8.0 keV). This flux is lower than 
the Suzaku observation by more than a factor of 100. 

2.2 Suzaku Observation 

Our observations of VII Zw 403 used the three functioning X- 
ray Imaging Spectrometer (XIS) cameras aboard Suzaku, one of 


2.3 Suzaku Timing Analysis 

Short term variability (e.g., quasi-periodic oscillations (QPOs), 
pulsed emission) has been observed in some ULXs (e.g. 
Strohmayer & Mushotzky 2003; Dewangan et al. 2006; Feng et al. 
2010). These features allow for identification of certain types of ob¬ 
jects (e.g., pulsars, eclipsing binaries, etc.) and can constrain size 
and mass of the emitting object in ways the X-ray spectra cannot. 

Observations of VII Zw 403 were taken in the normal timing 
mode of Suzaku which has a timing resolution of 8 seconds. Light 
curves were extracted from the same 3' circular regions from which 
the spectra were extracted. Energy filtering was applied such that 
only photons with energies in the 0.3 — 8.0 keV range were ac¬ 
cepted. Using the HEASARC tool xselect, we produced light 
curves with 200s time bins for sufficient counts. The timing anal¬ 
ysis that was done without background subtraction showed no sig¬ 
nificant variability beyond Poisson fluctuations for the XIS0,3 de¬ 
tectors. Analysis of the background-subtracted light curves still 
showed no significant variation in the XIS0,3 detectors. The XIS0,3 
detectors had net count rates of 0.25 and 0.26 counts/s, respec¬ 
tively. The sample and expected variances for these detectors were 
~ 1.4 X 10“^ and are consistent with a constant light curve. 

Assuming the count rate is constant, the XIS 1 detector exhib¬ 
ited variances that were much larger than expected from Poisson 
noise alone as indicated by an F-test of the variances with 1-P(F- 
test) > 0.99. The XISl detector is more sensitive to background 
variations than XIS0,3 because it is back-illuminated. The back¬ 
ground light curve for XIS 1 showed the same strong variation as 
the source light curve for XIS 1 and thus the background variations 
are the likely cause of the excess variance. To try and remove any 
effects due to background variation, we subtracted the background 
counts on a bin by bin basis for the source region. For a given time 
bin, we found the number of background counts in the background 
region, scaled this number by the ratio of the areas of the source 
and background regions, and subtracted this from the total counts 
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in the source region. The errors were added in quadrature with a 
weighting to account for the difference in extraction region area. 
This process removed some of the strong variations we were seeing 
from the XISl source region. The variations in the XISl detector 
were reduced but still showed excess variance with a signihcance 
of l-P(F-test) = 0.99. The net count rate for the XISl detector was 
found to be 0.28 counts/s with a sample variance of 2.0 x 10“® and 
an expected variance of 1.7 x 10“®. Due to the noise in the XISl 
detector, we exclude it from further timing analysis. 

Power spectra were calculated for each of the FI detector light 
curves using the method outlined by Leahy et al. (1983). Due to 
Suzaku’s orbit, there are periodic gaps in the data that occur as a 
result of either Earth occultation or passage through the South At¬ 
lantic Anomaly (SAA). These periodic gaps produce large peaks 
in the FFT of the light curves. To avoid this problem, we calculate 
power spectra for each continuous set of light curve data and com¬ 
bine the power spectra in the way described by Leahy et al. (1983). 
For each detector, 27 power spectra were averaged over a frequency 
range of 0.37 — 62.1 mHz (T = 2720 s) using 169 frequency 
bins. The power in the averaged power spectrum is distributed as 
a random variable with 2N degrees of freedom, where N is the 
number of power spectra that were averaged together. We found 
a maximum power of ~ 3.0 (Leahy) between the two detectors. 
The chance probability of detecting this maximum power in any 
frequency bin is p = 0.76 and is therefore not significant. 

2.4 Suzaku Spectral Analysis 

Due to uncertainty in the normalization between the BI and FI de¬ 
tectors, we allowed for a constant factor between the XISl and 
XIS0,3 spectra. This constant factor increases the overall normal¬ 
ization of the XIS0,3 spectra by about 5 per cent. 

We began by htting only single-component models that in¬ 
cluded interstellar absorption. We dehne an acceptable fit as hav¬ 
ing a value that lies in the lower 95 per cent of the respective 
X^-distribution. As summarized in Table 2, the single multicolored 
disk (DISKBB in XSPEC), simple power law models (POWERLAW 
in XSPEC), and the exponentially cutoff power law (CUTOFF in 
XSPEC) did not fit the data well. A broken power law or Comp- 
tonization model described the data within acceptable x^ values. 
The broken and cutoff power law models both provided signifi¬ 
cantly better fits than a single power law model. This is an indica¬ 
tion of spectral curvature, a ubiquitous feature in high quality ULX 
spectra. 

The Comptonization model provides the best fit of the single¬ 
component models and is the only model that can provide a direct, 
physical interpretation of the data. It exhibits an extremely cool in¬ 
put photon temperature of To < 0.09 keV with almost no intrinsic 
absorption and an electron coronal temperature of 2.15 ± 0.06 keV 
with a coronal optical depth of 11.8 ± 0.2. These parameters are 
consistent with other high-quality ULX spectra (Gladstone et al. 
2009). By looking at the residuals for these different fits, one can 
see that the vast improvement is due to the modelling of a clear 
break in the spectra around 5 keV (see Figure 2). The residuals also 
indicate that below about 0.8 keV there is a soft excess that cannot 
be fit by a single component. This provides the evidence for the 
second main feature of ULX spectra, a soft excess, and motivates 
fitting using more complex models. 

The appearance of a soft excess in the residuals of our single¬ 
component model fits motivates the addition of a second com¬ 
ponent to the models. The spectra we analyzed have sufficient 
net counts (> 10 000) such that multi-component models may 


be justified. We fit the spectra with a disk plus powerlaw model 
(diskbb+powerlaw), and we find two minima in the fits 
(Models 7 and 8). For models (7), (8), and (9) it is assumed that 
the disk is providing seed photons to the mechanism driving the 
power law. The fit for which the power law models the soft compo¬ 
nent and the disk models the hard component (Model 7) has been 
discussed in other papers (Foschini et al. 2004; Roberts et al. 2005; 
Roberts 2007) and was shown to be unphysical. For this reason, the 
diskbb-l-powerlaw model (Model 7 in Table 2) is not physically re¬ 
alizable. We also find that Model (8) does not provide an adequate 
fit due to the curvature above 5 keV and is improved by the use of 
an exponentially cutoff power law (Model 9). 

Gladstone et al. (2009, 2011) provide an observationally de¬ 
fined state for a set of high quality ULX spectra called the ultralu- 
minous state. This state is characterized by a soft excess modelled 
by a cool disc with a power law tail and a break in the spectrum 
above ~ 3 keV, similar to the break observed with Suzaku from 
VII Zw 403. In order to compare our spectra with the ULX sample 
of Gladstone et al. (2009), we use the diskpn model for the disk 
component as opposed to diskbb. The two disk models are con¬ 
sistent with each other to within 5 percent. Referring to Table 2, 
one finds that most of the two-component models provide ade¬ 
quate fits, but only the disk plus Comptonization model provides 
a physically realistic description of the spectra. In this model, the 
disk temperature is tied to the input soft photon temperature of the 
Comptonization model. The parameter values for this model agree 
with the compTT model without the disk component. However, the 
two-component model provides a fit that is significantly better (1- 
P(F-test) > 0.99) than the compTT model on its own. This is due to 
improved modelling of the soft excess by the additional disk com¬ 
ponent. 

From the Comptonization models, we find that the fits de¬ 
scribe a cool, optically thick corona (fcTe = 2.2 keV and r = 
11.4) consistent with values found by Gladstone et al. (2009) (i.e., 
kTs ~ 1 — 3 keV and t ~ 6 — 80). The very cool disc temperature. 
To < 0.1 keV, is consistent with disc temperatures predicted for 
intermediate mass black holes (To < 0.5 keV). However, this is 
unlikely to be representative of the actual temperature as the pres¬ 
ence of an optically thick corona would lead to obscuration of the 
accretion disc (Kubota & Done 2004), making direct observation 
of the inner accretion disc impossible. Any physical interpretation 
based on the inner accretion disc temperature parameter would be 
questionable and we simply use this parameter as a means to com¬ 
pare to the results of Gladstone et al. (2009). 

Sutton et al. (2013) divided the ultraluminous state as defined 
in Gladstone et al. (2009) into three spectral subtypes: broadened 
disc, hard ultraluminous, and soft ultraluminous. The broadened 
disc state represents a state in which the accretion rate is near the 
Eddington accretion rate. The hard and soft ultraluminous states 
are thought to be super-Eddington accretion states. As observed 
with Suzaku, VII Zw 403 has a kTm < 0.5 keV and F < 2 
(diskbb+powerlaw. Model 8 in Table 2) which puts it in the 
hard ultraluminous state. Using the criteria of Sutton et al. (2013), 
we find the spectral hardness to be F{1 — 10 keV)/T(0.3 — 
1.0 keV) = 3.80, which is the ratio of the unabsorbed fluxes in the 
indicated energy ranges. This places VII Zw 403 within the same 
region as the other hard ultraluminous objects in the hardness- 
luminosity diagram of Sutton et al. (2013). 
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Figure 2. Single Component Models: Suzaku spectra of VII Zw 403 (XISO: black, XISl: green, XIS3: red) rebinned to a minimum of 15f7 or a maximum of 15 
channels for visualization. Left: Model (2) power law fit (see Table 2). Right: Comptonization (compTT). The power law models cleai'ly show a break in the 
spectra above 5 keV and the need for a soft excess at lower energies (E < 0.8 keV). The Comptonization model provides an acceptable fit, but is significantly 
improved by the addition of a soft disk component (see Figure 3). 


Table 2. VII ZW 403 spectral fitting results: Flux (absorbed and unabsorbed Fx), luminosity and fitting are all in the 0.4 — 10.0 
keV energy band. An unabsorbed flux of 7.0 x 10~^^ cgs converts to a luminosity of 1.6 x 10“^^ erg s“^, with an assumed distance of 
4.34 Mpc. All models use the tbabs absorption model with Galactic absorption Nh = 3.91 x 10^^ cm~^, which was found using the 
Golden tool (http : / /cxc . harvard . edu/toolkit/colden . j sp), and the tbvarabs model for the intrinsic absorption such 
thatZ = O.O62Z0. 


Model 

X^/DoF 


Fx 

Nh 

kT/Ec 

r/r (*>) 




(10~^^ cgs) 

cgs) 

(lO^l cm-2) 

(keV) 



(1) Diskbb 

1749.8/1256 

6.7 

6.8 

< 0.009 

2.99±0.06 

_ 


(2) Power law 

1460.2/1256 

7.2 

lA 

< 0.06 

- 

L13±0.01 


(3) Cutoff PL 

1381.6/1255 

7.0 

7.2 

< 0.021 

17±3 

0.95±0.04 


(4) Broken PL 

1305.6/1254 

6.9 

7.1 

< 0.025 

5'3to.3 

L07±0.02, 1.7; 

+0.3 

-0.1 

Model 

X^/DoF 

pixhs (c) 

Fx 

Nh 

To 

kTe/Ec 

F/t 



(10~^^ cgs) 

(10“^^ cgs) 

(10^1 cm-2) 

(keV) 

(keV) 


(5) CompTT 

1263.8/1254 

6.8 

7.0 

< 0.1 

<0.09 

2.15±0.06 

1L8±0.2 

(6) Diskpn+compTT 

1251.1/1253 

6.9 

7.5 

L4±0.9 

0.09±0.01 

2.19±0.08 

1L4±0.4 



0.1 / 6.8 

0.4/7.1 





(7) Diskbb+PL 

1250.8/1254 

6.9 

7.2 

< 1.1 

^•^-0.1 

- 

3.0± 0.5 



6.4 / 0.5 

6.6/0.6 




(8) Diskbb+PL 

1404.2/1254 

7.2 

8.4 

2.7±0.7 



1.17± 0.02 



0.2/7.0 

1.0/7.4 




(9) Diskbb+Cutoff 

1258.9/1253 

6.9 

7.1 

< 0.14 

0.24±0.04 

6 4+^'^ 

0.55±0.! 



0.3 / 6.6 

0.3/6.8 





kT is the inner disk temperature for disk models. Ec is the break/cutoff energy for power law models, 
r is the photon index for power law models. For compTT, r is the plasma optical depth. 

For the two-component models, individual flux contributions are given below the total flux values using the convention: 

(disk component flux) / (PL or Comp. flux). 

To is the inner disk temperature. 

kTe is the plasma temperature for the compTT models and Ec is the cutoff energy for the exponentially cutoff power law models. 


3 DISCUSSION & CONCLUSIONS 
3.1 Spectral evolution of VII Zw 403 

Observations of VII Zw 403 have shown luminosities in the 1.3 — 
23 X 10^® erg range [0.3 — 8 keV]. Table 1 and Figure 1 show 


the known X-ray flux history of VII Zw 403. While in the low flux 
state, both Chandra and XMM-Newton observed the X-ray source 
to have a spectrum consistent with a hard power law (F = 1.75 — 
2.1). Our analysis of a recent Suzaku observation shows that the 
X-ray source transitioned to a flux more than 100 times brighter 
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Figure 3. Suzaku spectra of VII Zw 403 (XISO: black, XISl: green, XIS3: 
red) rebinned to a minimum of 15cr or a maximum of 15 channels for visu¬ 
alization. The model (black line) consists of an absorbed accretion disk (red, 
dashed line) and a Comptonization component (blue, dash-dotted line). The 
X-ray spectrum of VII Zw 403 is quite flat and curvature can be seen at 
energies > 5 keV. These spectral properties, along with the observed lu¬ 
minosity of 1.7 X 10"^® erg s“^, are consistent with ultraluminous states 
found in other extragalactic sources (Gladstone et al. 2009). 

(1.7 X erg s“^) than the XMM observation. More recent Swift 
observations indicate that the X-ray source has transitioned out of 
the high flux state. 

The high quality of the Suzaku spectra allowed for the iden¬ 
tification of spectral curvature above an energy of S.SIqj keV us¬ 
ing a broken power law fit (Model 4 in Table 2), a feature seen in 
many high quality ULX spectra. Using a disk plus Comptonization 
model, we fit the spectra and found our parameters to be consis¬ 
tent with the sample of ULXs in the ultraluminous state studied by 
Gladstone et al. (2009). This is consistent with the idea that a com¬ 
pact object with an accretion disk is producing photons that are 
being upscattered by an optically thick corona. 

Motch et al. (2014) have recently reported that the X-ray 
source P13, a ULX in NGC 7793, exhibits all the canonical prop¬ 
erties of ultraluminous sources and is in a binary system for which 
optical modulations constrain the black hole mass to be less than 
15 solar masses. This X-ray source has also been observed to 
vary in its X-ray luminosity by at least a factor of 40. Their fit¬ 
ting of a disk plus Comptonization model results in a coronal 
electron plasma temperature of = 2.11 q j keV and an op¬ 
tical depth of T = which are consistent with our re¬ 

sults for VII Zw 403. Other studies (e.g. Poutanen et al. 2007; 
Gladstone et al. 2009; Feng & Soria 2011; Kawashima et al. 2012; 
Sutton et al. 2013) also support the interpretation that similar ultra- 
luminous sources are StMBH binaries and not the more massive 
IMBH binaries. Drawing conclusions about VII Zw 403 from the 
results of Motch et al. (2014) and others would be speculative but 
a simple comparison may suggest that VII Zw 403 is a stellar mass 
black hole undergoing supercritical accretion. More uniform mon¬ 
itoring over a longer interval is needed to adequately constrain the 
luminosity evolution of the source. 


3.2 Alternative possibilities 

Within VII Zw 403 there is a second known X-ray source (X-2) 
that lies within the beam of Suzaku. This source has been observed 
to have a luminosity 25 times dimmer than the main source (X- 


1) (Brorby et al. 2014) we have discussed in this paper. These two 
objects are separated by a projected distance of 650 pc or about 
30 arcsec. The relative offsets of the Suzaku source with respect 
to the well-defined Chandra sources are 20 '.'8 for X-1 and 13'.'2 
for X-2. Suzaku’s absolute pointing has an uncertainty of about 
20'.' (Uchiyama 2008). The PSF of the XRTs has a broad half¬ 
power diameter (HPD) of ~ 2', but a sharp Gaussian core of width 
~ 10'.'. Given these two independent errors, the centroid position of 
the Suzaku source is consistent with both X-1 and X-2. It is there¬ 
fore possible that the dimmer source has brightened, but this would 
require a factor of 1000 change in luminosity, if due solely to this 
second source. The XMM and Chandra observations are able to re¬ 
solve these two sources and we find that the flux ratio between X-1 
and X-2 decreases by a factor of about 2.7 over the two observa¬ 
tions which is consistent with the factor of g change in the flux 
of X-1 alone (see Table 1 ). Thus, X-1 seems to be the source that is 
varying, whereas no significant variability is measured for X-2. 

From the study presented in this paper, there is no evidence 
that the main source in VII Zw 403 is a background AGN. The 
number of expected background sources within the D 25 ellipse 
coincident with VII Zw 403 is < 0.02, as determined by the 
log N — log S curves of Georgakakis et al. (2008). Thus, it is un¬ 
likely that either of the two sources associated with VII Zw 403 are 
background AGN. 

3.3 Comparison to I Zw 18 

Kaaret & Feng (2013) reported a transition from low to high flux 
for an X-ray binary in the lowest known metallicity BCD, I Zw 18. 
The X-ray source was observed to transition from a low-flux state 
with a hard power law spectrum to a high-flux state with a thermal 
or exponentially cutoff spectrum. Figure 3 of Kaaret & Feng (2013) 
shows that at high-flux the source exhibits curvature over the entire 
0.3 — 10.0 keV range. This type of high-flux state is consistent with 
Galactic black hole binaries and lies in stark contrast to the high- 
flux state of VII Zw 403. As one can see from Figure 3, VII Zw 403 
shows a relatively fiat spectrum and curvature only above 5 keV, 
which is more consistent with high-quality ULX spectra. 

The difference in the environments in which the X-ray sources 
of VII Zw 403 and I Zw 18 reside may provide an explanation 
for the differences in their spectral evolution. VII Zw 403 and 
I Zw 18 are dwarf galaxies with total stellar masses of 2 x 10^ Mq 
and 5.9 x 10® M© (Zhao et al. 2013), star formation rates of 

1.4 X 10“^ Mq yr“^ and 7.0 x 10~^ Mq yr“^ (Brorby et al. 
2014; Prestwich et al. 2013), and metallicities of Z/Zq = 0.062 
and Z/Zq = 0.019 (Izotov et al. 2007), respectively. In Figure 4, 
one sees that the source in I Zw 18 lies in a dense star-forming 
region, whereas the object in VII Zw 403 is far from the dynami¬ 
cal center, away from dense regions. Intermediate mass black holes 
are more likely to form in young, dense clusters where dynamical 
friction causes the runaway merger of stars over a relatively short 
timescale, tdf ^ 4 Myr (Portegies Zwart et al. 2004). The spec¬ 
tra of these objects (if in a binary system) are expected to behave 
as scaled-up versions of the Galactic black hole binaries. The fact 
that I Zw 18 exhibits this type of spectral state behavior and that 
it resides in a dense environment does not rule out the possibil¬ 
ity that it is an intermediate mass black hole binary. However, for 
the source in VII Zw 403, given its location and spectral shape, it 
seems unlikely that it could be an intermediate mass black hole. 
More observations that catch these objects in a broad range of lu¬ 
minosity states would be needed to conclusively determine their 
spectral state transition history and thus the nature of their primary 
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compact objects. Recent results by Bachetti et al. (2014) illuminate 
the need for more observations by showing that a known ULX in 
the nuclear region of M82, which contains many dense star clus¬ 
ters, is a super-Eddington accreting neutron star and not a massive 
hlack hole. 

Most extragalactic ULXs that exhibit spectral curvature are 
located away from the galactic nucleus and removed from very 
dense regions (Pakull et al. 2006; Roberts 2007; Grise et al. 2008). 
These ULXs, based on their location and spectra, are argued to 
be stellar mass black hole binaries emitting at super-Eddington 
luminosities, though the mechanism by which this happens is 
not completely understood. VII Zw 403 observations are more 
consistent with these properties, suggesting it may be a super- 
Eddington stellar mass black hole binary. The transient behavior 
of VII Zw 403 source is particularly interesting since most “stan¬ 
dard” ULXs are nontransient and may be active for up to decades 
in some cases (Kaaret & Feng 2009; Feng & Soria 2011). In this 
way, VII Zw 403 is more like the Galactic BHBs, spending most 
of its time in the quiescent state. More observations of the low and 
transition states may provide a means to understand the connection 
between the ultraluminous state and the states commonly seen in 
Galactic black hole X-ray binaries. 

3.4 Implications 

Study of nearby HMXBs can provide insight into the formation of 
high mass X-ray binaries formed in the early universe. Observa¬ 
tional and theoretical studies have shown that reduced metallicity 
enhances X-ray binary production and the total X-ray luminosity 
of a galaxy. The total luminosity, and hence the spectral shape of 
the emission from a population of HMXBs, is dominated by a rel¬ 
atively small population of sources in high luminosity states. To 
date, a total of eight low-metallicity BCDs have been found that 
contain ULXs (Thuan et al. 2014) as defined by having a luminos¬ 
ity > 10®® erg s“^ (Eddington luminosity of a 10 Mq black hole). 
Only two of these, VII Zw 403 and I Zw 18, have had sufficient 
counts to allow for a spectral study and produced inconsistent re¬ 
sults. Our discovery of this state transition of an X-ray source in a 
low-metallicity blue compact dwarf galaxy shows that the sources 
may undergo a transition to a spectral state with high-energy cutoffs 
at high luminosity, similar to those seen in other ULXs. 

Other low-metallicity, star-forming dwarf galaxies, such as 
Holmberg II (Goad et al. 2006) and Holmberg IX (Stobbart et al. 
2006; Gladstone et al. 2009; Walton etal. 2014), contain known 
ULXs that exhibit spectra consistent with the ultraluminous state, 
as defined by Gladstone et al. (2009). The metallicities of these 
galaxies are low (12 -f log(0/H) < 8.0), but do not meet the crite¬ 
ria for being XMPGs (12 -f log(0/H) < 7.7), mostly occupied by 
BCDs. Low metallicity galaxies, such as Ho II and Ho IX, may be 
useful for studies of ULX spectral evolution with metallicity. 

These results motivate further spectral studies of luminous 
HMXBs in nearby analogues of early universe galaxies. This will 
allow for better estimates of the spectral shape of HMXBs in the 
early universe, which affect the ionisation morphology during the 
Epoch of Reionization via thermal feedback, and the time at which 
reionisation was completed. 
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Figure 4. HST images with X-ray source (red cross) and dynamical centers (red X, (Lelli et al. 2012; Simpson et al. 2011)). I Zw 18 (right) has a luminous 
X-ray source within a dense star' forming region of the galaxy. VII Zw 403 (left) contains a X-ray source that is not near a dense star forming region and far 
from the dynamical center of the BCD. 
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